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Abstract 

The present article discusses about the effect of a Lee- Wick partner infested radia- 
tion phase of the early universe. As Lee- Wick partners can contribute negative energy 
density so it is always possible that at some early phase of the universe when the 
Lee- Wick partners were thermalized the total energy density of the universe became 
very small making the effective Hubble radius very big. This possibility gives rise to 
the probability of a bouncing universe. As will be shown in the article that a simple 
Lee- Wick radiation is not enough to produce a bounce. There can be two possibilities 
which can produce a bounce in the Lee- Wick radiation phase. One requires a cold dark 
matter candidate to trigger the bounce and the other possibility requires the bouncing 
temperature to be fine tuned such as all the Lee- Wick partners of the standard fields 
are not thermalized at the bounce temperature. Both the possibilities give rise to 
blue-tilted power spectrum of metric perturbations. More over the bouncing universe 
model can predict the lower limit of the masses of the Lee- Wick partners of chiral 
fermions and massless gauge bosons. The mass limit intrinsically depends upon the 
bounce temperature. 
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1 Introduction 



Recently the idea of constructing a unitary, gauge- invariant, Lorentz-invariant and divergence 
free theory of Quantum Electrodynamics by invoking unusual partners of the standard par- 
ticles in the Lagrangian, originally proposed by Lee and Wick [U Ej , has been generalized to 
construct a Lee- Wick theory of the standard model of particle physics [3]. This model, while 
still suffers from unsettled issues such as quantum instability of the vacuum due to gravita- 
tional interaction [3] , attracts a lot of attention in the literature for several phenomenological 
benefits. Namely, the mass of the Higgs field in such a quantum field theory is perturbatively 
stable against the quadratically divergent radiative corrections and thus the well-known "hi- 
erarchy puzzle" of the Standard Model of particle physics could be avoided. There are many 
other interesting applications of the Lee- Wick idea of which some are discussed here. In [S] 
a minimal extension of Lee- Wick Standard model (LWSM) is considered, in [6] a LWSM 
with more than one Lee- Wick (LW) partner for each standard model particle is studied. 
Ref. [7] deals with gauge-coupling unification in Lee- Wick framework and in [8j [9] analysis of 
two-Higgs doublet models where one of the doublet contains Lee- Wick fields is presented. In 
Ref. PU] the process gg — > h — >■ 77 is studied in the framework of LWSM where the authors 
predict small changes in the rate of these processes, due to the presences of Lee- Wick fields, 
from those rates calculated from other models such as universal extra dimensions. In [TTj 
Higgs pair production processes gg — > h h and gg — > h po are studied in LWSM framework. 

It was recently realized in [12] that a model constructed out of a Lee- Wick type scalar 
field theory is capable of yielding a bounce during the evolution of the universe within the 
framework of a homogeneous and isotropic cosmological background. In this model the 
energy density and pressure of all other kinds of fields are supposed to be negligible. In 
this scenario the Lee- Wick partners, which arise from higher derivative operators, evolve as 
a tracking solution to the normal matter (the standard scalar field) and can break certain 
energy conditions when the energy scale of the universe becomes high enough. Therefore, a 
non-singular bounce takes place which avoid the big bang singularity widely existing in the 
standard FRW cosmologies. Interests in such Lee- Wick non-singular bouncing models have 
increased recently as these models [131 EH] can be considered as an alternative to inflation and 
can be used to explain the origin of large scale structure in the universe. In this paradigm 
the primordial vacuum fluctuations leaves the Hubble radius during a matter-dominated 
contracting phase and then form a nearly scale-invariant power spectrum after the bounce. 

However, it was found that non-singular Lee- Wick bouncing cosmologies induced by 
pressure-less matter are in general unstable while accompanied with radiation and a bouncing 
solution might be achieved only when an extremely fine tuning of the initial phases of the 
field configuration is assumed [H)J US]- This problem can be roughly understood as follows. 
The normal matter which scales function of the scale factor a could be easily 
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suppressed by radiation matter which scales as a -4 in the contracting phase. If we simply add 
a radiation component to non-relativistic matter within a framework of a Lee- Wick bounce 
model without fine tuning of the initial values of radiation then the radiation would become 
dominant, during the contraction phase of the universe, and produce a big crunch singularity 
before the bounce takes place. In Ref. [15], the effects of Lee- Wick partners of radiation 
fields were taken into account as well but the energy scale of the universe is assumed to be 
below the scale of thermalization of Lee- Wick radiation. Then, it was found the effective 
energy density of Lee- Wick radiation increases slower than that of regular radiation and thus 
cannot help to realize the bounce at low temperature regime. 

But unlike these models of bouncing cosmologies, where the energy density and pressure 
of the universe arises from a scalar field (whose energy scales as a -3 ) or a radiation field 
(whose energy scales as a -4 ), one can conceive a phase of the early universe where we have 
the normal standard model particles and their Lee- Wick parters (as virtual resonances) in 
thermal equilibrium as shown in Refs. [T71 [18]. This phase is a non-standard radiation 
dominated phase where one can have a state parameter w ~ 1 and consequently the energy 
density scales ~ a -6 . It is explicitly described in Ref. [18] that in such kind of a Lee- Wick 
partner infested radiation dominated universe the time-temperature relation in standard 
cosmology changes and the universe may go through out-of equilibrium processes when 
some of the Lee- Wick partners become non-relativistic. It was noticed that these kind 
of cosmologies can produce negative energy densities and pressure when one includes chiral 
fermions and massless gauge bosons [T9|[20]. This particular point turned out to be a serious 
drawback for these models. It was shown that only under some specific conditions these kind 
of an early radiation domination can exist. On the other hand bouncing cosmologies naturally 
require a time (bounce time) when the Hubble parameter turns out to be zero (indicating a 
zero energy density). It turns out that one can construct an interesting model of a bouncing 
cosmology from the non-standard Lee- Wick radiation phase. 

In this paper we study the possibility of realizing a Lee- Wick bounce at high temper- 
ature by taking into account the radiation fields. The Lee- Wick resonances are produced 
dramatically when the temperature of the radiation fluid is much higher than the mass of 
heavy Lee- Wick fields. As the contribution of Lee- Wick radiation to the total energy den- 
sity is always negative it enlarges the viable phase space for the universe to experience the 
nonsingular bounce. Unlike the method developed in [J5] in which the radiation fields are 
described by the field configuration of gauge fields, we directly study thermodynamics of 
normal particles and their Lee- Wick partners following [T7[ US] . However, for a radiation 
dominated universe purely dominated by standard fields and their Lee- Wick partners, one 
still cannot obtain a nonsingular bounce. This is because, at the bouncing moment when 
the energy density of Lee- Wick resonances are able to cancel that of normal radiation, the 
pressure of the universe keeps positive and then prevent the background to evolve into an 
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expanding phase. Thus we introduce a cold dark matter (CDM) component in our model 
which can easily solve this problem. As a consequence, our cosmological model, which real- 
izes a nonsingular bouncing solution in the thermal Lee- Wick theory, consists of cold dark 
matter fluid, radiation, and Lee- Wick resonances. 

If one does not include the CDM component then one cannot obtain the required condi- 
tion of the bounce as discussed in the last paragraph. This fact is true if all the standard 
particles and their Lee- Wick partners are in thermal equilibrium near the bounce time. But 
if it happens that the bounce temperature is such that some of the Lee- Wick partner's mass 
are more than the bounce temperature then the plasma around the bounce time will have 
less Lee- Wick partner. Interestingly, this kind of a model where the bounce temperature is 
less than some of the masses of the Lee- Wick partners can give rise to a non-singular bounce. 
This model does not require any CDM candidate to make the pressure negative at the bounce 
time. In the present article we will also discuss about this alternative bouncing model. Both 
of the bouncing models yields a blue-tilted power spectrum of metric perturbations in the 
expanding phase after bounce. 

The outline of this paper is as follows: In Section 2 we briefly introduce the thermody- 
namics of Lee- Wick particles. In Section 3, we study the condition for a Lee- Wick thermal 
bounce and then provide a concrete cosmological model which involves both kind of bounce 
scenarios as described above. Section 4 is devoted to the discussion of cosmological pertur- 
bations in both of these models. The final section presents some conclusions and discussion. 
For completeness we have attached an appendix at the end of the article which discusses 
about the mode matching conditions near the bounce point. 

2 Thermodynamics of Lee- Wick particles revisited 

In this section we will give a brief account of the thermodynamics of Lee- Wick particles. We 
will show that in a realistic model of Lee- Wick theory, where all the degrees of freedom of 
the standard particles along with their Lee- Wick partners are taken into account properly, 
the energy density of the Lee- Wick resonance dominated universe would become very small 
when the heaviest of the Lee- Wick partners thermalize due to the very high temperature 
of the universe. This property is in favor of a bouncing scenario which requires the Hubble 
parameter (which is directly proportional to the square-root of the total energy density of 
the universe) to vanish at the time of bounce. 

It is noted in [T7J [18] that the Lee- Wick particles behave very differently from their 
standard partners in the relativistic regime. For the standard bosons and fermions the 
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energy density, pressure and the entropy density read as 



(sm) _ ff7r 2 T 4 (sm) _ ^tt 2 T 4 (sm) _ 2gn 2 T 3 

Pb ~ 30 ' Pb ~ 90 ' Sb ~ 45 ' 1 j 

and 

{sm) _ 7gn 2 T 4 (sm) _ 7g7r 2 T 4 (sm) _ 7grc 2 T 3 

Pf ~ 2 4o ' p f ~ 720 ' S / " 180 ' 1 j 

respectively. Here the subscripts b and / stands for bosons and fermions and g stands for 
any number of internal degree of freedom of the relativistic species. On the other hand, it is 
shown in [T7J [18] that the Lee- Wick particles at high energies contribute negatively to the 
energy density, pressure and entropy density as 

(lw) M 2 T 2 \ 

Pb = -g[— , (3) 

(lw) M 2 T 2 \ 



for the bosonic Lee- Wick partners and 

pf W) = -9 ( - ) > ( 6 ) 

H LW) = ( ^5n- - ) > ( 7 ) 

(LW) 

s) = -9 

for the fermionic Lee- Wick partners. In the above equations M is the mass of a generic 
Lee- Wick partner and as the system is relativistic T 3> M. As the Lee- Wick particles were 
initially introduced as resonance particles to overcome the divergences appearing in a theory 
and consequently when T ^> M the standard model particles are naturally ultrarelativistic. 
If one considers a toy model where each standard particle is accompanied by one Lee- Wick 
partner and the number of degrees of freedom of the Lee- Wick partner is same as that of 
its standard partner then in such a scenario the total energy density, pressure and entropy 
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density turn out to be positive as [T7J HB] 

Pb 
Pb 

for the bosonic sector and 

Pf 
PS 

for the fermionic sector. But this simple scenario, where each standard particle is accompa- 
nied by one Lee- Wick particle with equal number of degrees of freedom, does not serve the 
purpose while dealing with realistic scenarios. It is discussed in [19] that in the fermionic 
sector each chiral fermions require two Lee- Wick partners to eliminate the higher derivative 
terms in a initial higher derivative Lagrangian to construct a Lee- Wick theory. Such dis- 
crepancies also arise in bosonic sector when one considers massive Lee- Wick partners (with 
3 degrees of freedom) of massless gauge bosons (with two longitudinal degrees of freedom) 
[20] . Both these cases lead to unacceptable scenarios with negative energy density. It is 
shown in [20] that considering both the fermionic and bosonic sector together one still can 
come up with a scenario where the total energy density, pressure and entropy density can 
become positive. 

In a realistic scenario, where each fermion is accompanied by two Lee- Wick partners and 
the extra degrees of freedom of the massive Lee- Wick partners of each massless gauge boson 
are taken into account, the total energy density of the realistic Lee- Wick plasma can be 
written as [20] 

P = -24 9 * nT - 240^ T - 30 nT 
Here the new number of degrees of freedom g^^ is given as 



g*N — 
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(sm) . (lw) gM 2 T 2 
»f + »f --48- 

(sm) (LW) _ gM 2 T 
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(12) 
(13) 
(14) 



(15) 



where for bosonic particles stands for the number of internal degrees of freedom ^ 
for the partners of massive standard bosons (may be 2 or 1), while for standard massless 
vector boson partners it equals gi + 1 where primarily g± — 2. Also the unpaired fermionic 
contribution comes with gp where 

9f= Yl 9iF ( ~£ I ' ( 17 ) 



i=fermions 



where gp solely arises from the unpaired fermionic Lee- Wick partners of the standard model 
particles. The quantity n is defined as 



n = 

«=massivc vcct. bosons 



T 



Here Tj is the temperature at which the 2 th massive Lee- Wick vector boson partner of a 
standard massless gauge boson is equilibrated and n denotes the number of massive vector 
boson partners of massless standard gauge bosons if all the species are in thermal equilibrium 
at the same temperature T. The sum appearing in Eq. (TT8]) does not include all the massive 
Lee- Wick vector boson partners, but includes only those which are partners of massless 
standard gauge bosons. The total pressure of these Lee- Wick infested plasma is given by 

p= Q*n-l Qf-l nl . 19 

v 2A y*N 720 90 V ; 

It must be noted that if the energy density of the Lee- Wick partner infested universe is 
positive then the pressure of the same universe must be positive. It is also to be noted 
that the negative contributions to the total energy density and pressure are coming from the 
Lee- Wick partners of the chiral fermions and massless gauge bosons. 



3 Model(s) of Lee- Wick thermal bounce 

Our main aim is to construct a non-singular bouncing cosmology within the arena of thermal 
Lee- Wick scenario in which the universe initially starts its evolution in a contracting phase 
and evolves into a standard thermal expanding phase smoothly and continuously through 
a non-singular bounce. We first start with a general discussion on the conditions required 
for a non-singular bounce and would consider a simple model where the universe consists 
of bosonic and fermionic radiation fields and their thermalized Lee- Wick resonances. We 
will show that this simple toy model, which consists of only radiation fields accompanied 
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with their thermal Lee- Wick partners, is incapable of achieving the conditions required for 
a non-singular bounce and then we will try to illustrate two possible modifications to this 
simple model where a bounce can be achieved. 
We start with a spatially flat FRW metric 

ds 2 = dt 2 - a 2 {t)dx 2 , (20) 

where a(t) stands for the scale factor. The dynamics of the FRW universe is described by the 
Hubble parameter H = a/a and its time derivative H which obey the well-known Friedman 
equations: 

H 2 = ^P, (21) 
H = -47rG(p + p) . (22) 

At the moment of the non-singular bounce, i.e. at t = t B , the Hubble parameter vanishes 
Hits) — and in order to ensure that the universe enters an expanding phase one also 
requires Hits) > at the same moment. 

Let us consider that in the contracting phase the universe was initially radiation dom- 
inated with the standard particles (with w — 1/3). As the scale factor contracts the tem- 
perature of the Universe increases and the Lee- Wick particles (which are heavier than their 
standard model partners) will gradually start to contribute thermally as t — > CT and just 
before the bounce the energy density of the universe would be given by Eq. ( fl5l) . To satisfy 
the first condition of bounce i.e. H(t B ) = one requires p(t B ) = which can be seen from 
Eq. (liZTj) . Thus, assuming that just before the bounce all the particle species and their Lee- 
Wick partners are thermalized, it can be seen from Eq. ([15]) that at the time of bounce the 
energy density vanishes i.e. pits) = to yield 

M 2 ~ ^2 77T 2 4 7T 2 . 

^9*nT b = —g F T B + —nT B , (23) 

where T B is the temperature of the universe at bouncing time t B . This shows that at bounce 
the pressure given in Eq. (fl9|) would be positive as 

Pits) = t^~9fT b + n -nT B > . (24) 
Thus at the bouncing point one has from Eq. (122]) that 

H(t B ) = -47iG(p(t B ) + p{t B )) < , (25) 
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which indicates that after the bounce the universe fails to enter into an expanding phase. 
This means that in this simple cosmological model which is comprised of pure radiation 
and their thermalized Lee- Wick resonances, one cannot obtain a bouncing solution for a 
universe. Now we will discuss two distinct cases which slightly deviate from this simple 
scenario discussed above and where a bouncing universe scenario can be realized. 

3.1 Lee- Wick infested radiation plasma and a CDM component 

First we consider a CDM component along with the radiation plasma which is composed of 
thermal standard particles and their thermalized Lee- Wick partners. The CDM candidate 
is non-relativistic throughout the whole evolution and hence we can neglect its pressure in 
the following analysis. The CDM component too can have Lee- Wick partners according to 
the Lee- Wick scenario, but these partners would be heavier than the CDM particle under 
consideration. Hence the Lee- Wick partners of the CDM component could never thermalize 
in the evolution (as the CDM remain non-relativistic throughout) and would remain off-shell. 
Thus these Lee- Wick partners of CDM component would not contribute thermally to the 
energy density of the radiation plasma. Thus the energy density and pressure of the CDM 
component would be 

p D = n D m , p D ~ , (26) 

respectively, where hd and m are the number density and mass of the CDM particle respec- 
tively and m ^> T during all the phases of evolution of the universe. 

Let us briefly illustrate the scenario under consideration. We naively imagine an inverse 
picture of evolution of the expanding phase of the universe in the contracting phase. Hence in 
the contracting phase the universe initially evolves in a CDM dominated phase with w — 0. 
As the universe contracts the temperature of the universe increases and the universe enters a 
standard radiation dominated phase with w — 1/3. As the temperature of the universe keeps 
on increasing with the contraction the Lee- Wick resonances start to contribute thermally and 
negatively to the energy density of the radiation bath. In that case, though the temperature 
of the universe increases with contraction the energy density starts to decrease with more and 
more Lee- Wick particles (mainly the partners of chiral fermions and massless gauge bosons) 
contributing thermally. The CDM component, during this period, evolve non-relativistically 
with the background evolution decoupled from the thermal bath. Once the energy density 
of the Lee- Wick infested radiation bath becomes slightly smaller than the energy density of 
the CDM component, the later dominates the energy density of the universe and helps a 
bounce to occur, as will be shown next. This is not a conventional matter domination phase 
as during this period the state parameter of the cosmic fluid would not be zero. This whole 
scenario is shown graphically in Fig. [TJ 
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Figure 1: Figure showing the bounce point and the energy balance of the universe. The 
CDM domination near the bounce point is due to the addition of the non-standard CDM 
like component in the theory. 



According to the above discussion, the total energy density of the universe in this scenario, 
which consists of the energy density of the Lee- Wick infested radiation plasma pii(t B ) and 
that of the CDM particles Pd(^b), at the bouncing time t B turns out to be 

p(t B ) = PR.{t B ) + PdQb) = ^9*nTb - 7^9fT b - 7^nT B + p D (t B ) , (27) 

whereas the expression for the pressure of this Lee- Wick infested radiation plasma aided with 
the CDM particles remains the same as given in Eq. ffl9|) . If one wants to make p(t B ) < 0, 
required for a bounce at a temperature Tb, then the first constraint on the model is given 
by 

M\ 2 77T 2 _ 4 7T 2 4 

"24 9 * nT * < 720^ Tb + 90" Tb • (28) 

If this constraint on the system is fulfilled, then one can indeed have a bouncing universe 
where the energy density of the CDM candidate has to satisfy the second constraint: 

PD(t B ) > T^gpT^ + ^nT^ (29) 
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which is obtained by requiring that p(tg) = and using the first constraint given in Eq. 
Along with these two constraints we also must have the other constraint about the relativistic 
nature of the heaviest LW particle at Tg which reads as 

T B > M, (30) 

This condition along with Eq. yield the constraint on g* N as 

7tt 2 4tt 2 , , 

9*n > —9f + —n. (31) 
30 15 

If we fix the free parameters of this model as Tg ~ 10 16 GeV and roughly use the 
standard model values for calculating the number of internal degrees of freedom as, <jp ~ 40 
and n ~ 10, then these three parameters yield Pd(^b) > 4 x 10 64 GeV 4 and g*N > 100. Also 
to get pits) = one requires from Eq. (|27|) Pr^b) = — Pd^b) which yields M ps 9 x 10 15 
GeV which is the mass of the heaviest Lee- Wick particle in the model under consideration. 



3.2 Bouncing Universe within thermal LW framework without a 
CDM component 

We can have another scenario where only radiation can trigger the bounce. But for this to 
happen some conditions have to be met. The conditions are related to the mass of some 
Lee- Wick partners of the normal particles and the bouncing temperature Tg. In this case 
we consider that the bounce temperature Tb is smaller than the masses of some Lee- Wick 
particles. If this happens then some standard particles' Lee- Wick partners will not be able 
to thermalize as their masses are higher than Tg, but their standard partners being lighter 
would be thermalized at Tg. Due to this effect these (unpaired) standard particles will 
contribute only with positive radiation energy density. Such standard particles whose Lee- 
Wick partners could not thermalize at Tg can both be bosonic and fermionic. In this scenario 
the total energy density of the Universe can be written as 

M 2 7?r 2 tt 2 r 

Pits) = -^~9*nT 2 b - l —g F T B - -nT B + |t 4 , (32) 

and the total pressure would be 

p(*b) = ^9*nT 2 b - t^9fT b - ^nT 4 + -T% . (33) 

Here c and d encapsulate all the number of degrees of freedom and the corresponding factors 
of those standard particles whose Lee- Wick partners could not thermalize at bouncing point 
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and M is the mass of the heaviest Lee- Wick particle which could thermalize at or before 
Thus in this case too we have M < Tb- To satisfy the bouncing condition p{ts) < one 
requires (using Eq. ( 133]) ) 

M\ 2 77T 2 _ 4 7T 2 4 L 

Q*n1r< Qf1r-\ nin Jr, 34 

24 B 720 B 90 B 4 B K J 

and to achieve the other bouncing condition i.e. pits) = one gets 

c — d > - — Gf^r H nTt , (35) 

90 y B 45 B y ' 

where we have used Eq. (|32|) and Eq. (134|) . This scheme of a bouncing universe is also a 
possibility as for any thermalized bosonic or fermionic standard model particle p = \p and 
thus d = |c. Hence the above condition can be satisfied by choosing an appropriate set of 
parameters. 



4 Cosmological perturbations in bouncing universe 

In the previous section two distinct thermal Lee- Wick scenarios are considered where a non- 
singular bouncing universe scenario can be achieved. We devote this section to study the 
dynamics of linear cosmological perturbations generated in these models of thermal Lee- 
Wick bounce and to investigate the nature of power spectrum generated in such models. We 
focus on adiabatic fluctuations and consider matter components without anisotropic stress 
(we refer to [2T] for a comprehensive review on the theory of cosmological perturbations). In 
general, there are mainly two methods of analyzing cosmological perturbations in bouncing 
scenarios. One is to introduce a canonical variable of perturbation mode v (known as the 
Mukhanov-Sasaki variable) of which the quadratic action is of canonical form in the frame of 
conformal time coordinate. This variable is associated with the curvature fluctuation in co- 
moving gauge ( through v = z( where z is a background dependent coefficient and is roughly 
proportional to the scale factor if the background evolution is stable. A detailed calculation 
using a generalized canonical variable v in non-singular bounce cosmologies was carried out 
in [22] and it was shown that such a method is better suitable for those early universe models 
in which the primordial perturbations are originated from quantum fluctuations. Another 
method is to study the gravitational potential $ directly. The advantage of this method over 
the other is that it is much easier to impose initial conditions while analyzing gravitational 
potential $ as cosmological perturbations [23]. We will follow the second method for our 
analysis of cosmological perturbations. 
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Figure 2: Figure showing the contracting and expanding phases of the universe in the bounc- 
ing model. The horizontal axis is a comoving spatial coordinate and the vertical axis is 
conformal time. The main plot is of the Hubble radius \H\~ l and the wavelength A of fluc- 
tuations with comoving wave number k. kjjy is roughly the wave number at the onset of 
the bouncing phase whose value is dictated by the microphysics of the bounce. The heavy 
dashed curves on the main figure corresponds to the bounce phase. 

Before going into the details of the analysis of cosmological perturbations and the corre- 
sponding power spectra let us discuss briefly how the modes of these primordial perturbations 
evolve in the contracting and the expanding phases. For the discussion let us only consider 
the case with CDM component and Fig. [2] explains the essence of this particular bouncing 
model. Similar analysis also holds for the other case which does not consider any CDM 
component. In Fig. [2] it is shown how the Hubble radius \H\~ l behaves in the space-time 
diagram. The lower most heavy curve shows that the universe is undergoing a contracting 
phase where the energy density of the universe is changing from the standard radiation dom- 
inated phase (p oc T 4 ) to the mildly Lee-Wick radiation phase (p oc T 2 ). In this phase, the 
dark matter energy density is not important as it is mainly a radiation dominated phase. 
At the beginning of this Lee- Wick infested radiation phase, as all the Lee- Wick partners are 
not thermalized, the effective energy density will look like 

p ~ aT A + bT 2 - cT A , (36) 

where a, b and c are some constants whose exact value is not of much importance in this 
discussion. The first term aT 4 is the dominant term in the standard radiation domination 
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as shown in Fig. [T] and arises because the radiation domination is mainly dictated by the 
standard fields. The bT 2 term designates that some of the Lee- Wick partners have entered 
the scene and they have paired up with their normal partners to give the T 2 nature of energy 
density. The last term cT 4 is the contribution from the unpaired Lee- Wick partners of chiral 
fermions and massless gauge bosons. The schematic form of energy density, in Eq. (136]) . is 
similar to the form of energy density in Eq. (I15p except that here one also takes into account 
the standard thermal radiation part (glT 4 ). In presenting the energy density of Lee- Wick 
thermal radiation in Eq. (fl5|) one assumed that all of the standard field's Lee- Wick partners 
are in thermal equilibrium. In the initial stage of the contracting phase of the universe 
when the temperature of the system is not very high it is natural that most of the Lee- Wick 
partners would not have thermalized and so Eq. (136 p holds true. In this phase b 3> c and 
a c, which amounts to an assumption that most of the Lee- Wick partners of the chiral 
fermions and the massless gauge bosons (which give rise to the negative radiation energy 
density) are predominantly heavier than the other Lee- Wick partners. If this was not the 
case then the Lee- Wick partners of the chiral fermions and massless gauge bosons should 
have dominated the energy-density at the onset of Lee- Wick infested radiation era and the 
total energy density of the universe should have decreased dramatically. 

At the start of the contracting phase, when a > b, the state parameter w ~ 1/3 designat- 
ing a standard radiation phase. After some time b > a (and b ^> c) and the state parameter 
becomes w ~ 1. This phase is shown as the Lee- Wick radiation domination phase in Fig. 
[TJ As because of the b and c terms in the energy density the effective energy density of this 
phase rises but its rise is slowed down (as a result the growth of the Hubble parameter is 
inhibited) and consequently the Hubble radius goes on decreasing but with a slower rate. 
During this phase the modes with wavelength X = 1/k leaves the Hubble radius as shown in 
the Fig. [3 

As soon the conformal time is about the c term in the effective energy density starts 
to dominate and now c > a , b and the radiation part of the effective energy density starts to 
decrease rapidly. In this phase the radiation energy density starts to decrease as more and 
more Lee- Wick partners of the chiral fermions and massless gauge bosons thermalize and 
the Hubble parameter effectively decreases rapidly and the CDM energy density starts to 
dominate. This phase of the universe is represented by the non-standard CDM domination 
phase in Fig. [TJ The CDM contribution to the energy density ultimately balances the 
negative energy density due to Lee- Wick radiation and leads to a bounce. 

If we assume symmetry in time then a symmetric analysis can be presented for the 
expanding phase of the universe. The heavy dashed curves from to corresponds to the 
bouncing phase of the universe where the microphysics is mainly governed by the Lee- Wick 
partners of chiral fermions and massless gauge bosons and the dark matter sector. The solid 
heavy curves stand for the contracting and expanding phases of the universe. The modes 
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with wavelength A which went out during the contracting phase starts to re-enter the Hubble 
radius after in the expanding phase. In Fig. [2} tb is the time of bounce where as r% 
stands for the time where there should have been a Big Bang or Big Crunch like singularities 
in absence of the bounce. In Fig. [2J fcuv is the maximum possible wave number, which 
becomes super-Hubble, in such a bouncing scenario as it corresponds to the minimum radius 
of the possible Hubble radius possible. In absence of the bounce k\jy could have been much 
bigger, effectively of the Planckian order. 



4.1 Cosmological perturbations in contracting and expanding phases 
and their matching conditions 

Here we will discuss how cosmological perturbations generate in a contracting phase and 
evolve into the expanding phase in a generic bouncing universe scenario. The linearly per- 
turbed FRW metric in longitudinal gauge can be written as [21] 



ds 2 = a(r) 2 [(1 + 2$)rfr 2 - (1 - 2<$>)dx 2 ] , (37) 

where $ is the gravitational potential which characterizes metric fluctuations and r is the 
conformal time. At linear order the scalar metric fluctuations evolve independently. Thus 
we are able to study the evolution of $ by following one of its Fourier mode with a fixed 
comoving wave number k and the perturbation equation of $fc is given by 

+ 2 V H& k + (c 2 s k 2 - 2eK 2 + 2 V U 2 ) $ fc = , (38) 

where "H = a' /a and the prime denotes the derivative with respective to r. The sound speed 
parameter c s in the above equation is usually determined by thermodynamical property of 
the background system. Moreover, we have defined two background dependent parameters 
e = —H/H 2 and rj = —H/2HH. For a constant background equation of state (say w) these 
two parameters are equal and are totally determined by the background equation of state as 

e = v=l(l + w). (39) 

The situation which we are interested in here is a contracting phase of the universe 
dominated by radiation fields along with their thermalized Lee- Wick partners evolving in a 
stable background which evolves with a constant state parameter w. In that case the scale 
factor and the Hubble parameter of the Universe can be written as 

0~(T-7j,)^ , K~ I p- (40) 

(l + 3w)(r-r l B ) 
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respectively, where r B is some moment at which the Big Bang or the Big Crunch singularity 
would occur in absence of the non-singular bouncing phase. In our particular model t b 
t b ' t b as shown in Fig. [21 As a consequence, the perturbation equation given in Eq. 
can be simplified using Eq. fl39|) as 

^ + -^-& k + c 2 s k 2 $ k ^0 , (41) 

r - T B 

with 

1 + 3w 

Generically, the solutions to Eq. [41] are composed of two linearly independent Bessel func- 
tions. On super-Hubble scales they correspond to a constant mode (called as the D mode) 
and a time-evolving mode (called as the S mode). Specifically, the general solution to the 
perturbation equation on super-Hubble scales can be expressed as 

^ = D ± + S ± (^^-J , (43) 

where D and S are the mode coefficients and the subscripts ± represent the expanding 
and the contracting phase of the universe respectively. In the above expression tb is the 
conformal time at the bouncing point. It is also assumed in the above solutions that the 
background dynamics of the universe both in the contracting and in the expanding phases 
are governed by the same steady equation of state w. 

Now one needs to know how to transfer the primordial fluctuations generated in the 
contracting phase of the universe through the bounce. This issue was initially studied by 
replacing the bounce with a matching surface across which the perturbation modes are con- 
nected by using the Hwang- Vishniac [24] (or Deruelle-Mukhanov [25]) matching conditions. 
Later, it was found that for a non-singular bounce one can evolve the fluctuations through 
bounce both numerically and analytically [T2~| |2T)| |2"T] . Thus, relations between the mode 
coefficients in the expanding phase and those in the contracting phase can be calculated 
explicitly [12] (also see [23] for a general discussion^) . A general transfer relation between 



1 It is interesting to observe that for a nonsingular bounce model based on a closed space geometry, the 
primordial perturbation would be dramatically affected by the spatial curvature during the bounce even 
at large scale limit 28 . However, since in the model under consider the bouncing solution is achieved by 
introducing matter components which violate null energy condition in the frame of flat spatial coordinates, 
the curvature term will not be involved in the following transfer relation, as has been verified explicitly in 
the cold Lee Wick bounce model H2l. 



1(3 



these coefficients in the contracting and expanding phase can be written as 



D+ = 0{l)D- + 0(l)(-£-) SL, 



V fc uv/ 



(44) 



where fcuv is a normalization scale which is set by the micro-physics of the bounce. It 
corresponds to the inverse length scales near around in our model as shown in Fig. |2j A 
brief derivation of the above relation is presented in Appendix |A] 

4.2 Thermal fluctuations in the case which includes CDM 

In the model which we consider here the universe experiences a Lee- Wick infested radiation 
dominated period in contracting phase just before the CDM component takes over to yield 
a bounce. In the following we will study the generation of primordial power spectrum 
of curvature perturbation arisen from thermal fluctuations of Lee- Wick infested radiation 
plasma. Due to the existence of Lee- Wick partners, the heat capacity of the radiation 
plasma is different from the conventional one and thus affect the scale dependence of the 
power spectrum. Then we will study the stability of this cosmological system and estimate 
the mass bounds of the Lee- Wick partners. 

4.2.1 Generation of primordial power spectrum of thermal fluctuations in IR 



To determine the IR regime of the primordial power spectrum of the thermal fluctuations 
generated in the contracting phase of the universe, one should consider the modes which 
leave the horizon at earliest times of the contracting phase. From Fig. [1] we see that at the 
onset of radiation dominated era in the contracting phase the radiation fluid is dominated 
mainly by standard particles (as Lee- Wick partners being heavy thermalize at later stage) 
yielding w r ~ 1/3. The radiation energy density during that time thus can be written as 



where g sm takes into account the number of degrees of freedom of the standard particles 
thermalized initially. We also consider that the background during these initial phases of 
radiation domination still evolves under the influence of the CDM yielding w ~ and the 
index v — 5 introduced in Eq. (J42]). 

The correlation function of the energy density of such a cosmic fluid in thermal equilib- 
rium can be written as 



regime 



p r ~ g sra T 4 /4:, 



(45) 




(46) 



17 



where V is a fixed volume determined by the correlation length R of the thermal system 
which in a cosmological setup is roughly of the same order of the Hubble radius R ~ 1/H 
and a Fourier transform of any generic quantity <j) with the corresponding Fourier mode <p k 
is defined as 

^ x) = ^/(^ e!k '^ W - (47) 

The parameter Cy is the heat capacity of the radiation plasma and is defined in terms 
of the expectation value of the internal energy within the thermally correlated volume. 
These fluctuations in the energy density of the radiation fluid are coupled to the metric 
perturbations during the contracting phase through the time-time component of the 
perturbed Einstein equation which is given as 



- 3H(H$_ + $'_) + V 2 $_ = a 2 5p/2M 2 , (48) 

where Mp is the reduced Planck mass. When these perturbations leave the Hubble radius 
during the contracting phase (i.e. when the corresponding wavenumber of the mode is of 
the order of the comoving Hubble radius k ~ aH), the three terms contribute equally to the 
amplitude of $_ yielding 

, , a(r k ) 2 5p k ^ na\r k )C\l\r k )T{r k ) 

1 k 1 ~ 2k 2 M% ~ 2 1 /2fc7/2 M 2^( Tfc ) ' I 

up to a constant of order 0(1). We have used Eq. ( )48|) to derive the second expression and 
Tfc is the conformal time at the time of Hubble-crossing of the mode k. All the quantities in 
the above equation is to be derived at the time of Hubble crossing of the mode k. 

Now, the key issue is to find out the explicit form of heat capacity of the thermal system 
under consideration. By definition, the heat capacity of a thermal fluid is determined by 

C v = v %- ( 5 °) 

As has been mentioned before, the volume V is determined by the thermal correlation length 
which is in order of the Hubble radius. Consequently, using Eq. (145|) the heat capacity at 
the Hubble-crossing moment of the IR modes is given by 

Cy(r k ) « Y^ m (r fe )^P • (51) 

We note here that the perturbations in the radiation energy density scales with the scale 
factor w r as 

5p ~ a~ 3il+Wr \ (52) 
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where the entropy density of a fluid is defined as 

5p + 5p 5p 
s= — f — = {l+w r )—. (53) 

For an adiabatic expansion of the Universe the comoving entropy remains conserved which 
means a 3 s is a conserved quantity. Thus the temperature of the radiation plasma evolves as 

T ~ a~ 3Wr , (54) 

which follows from the previous two equations. Moreover, the scale factor of the universe 
evolves as a ~ r?'"" 1 ' in the contracting background which can be seen from Eq. (T4T)]) and at 
the time of Hubble crossing of the mode k one has Tk = — k- Then we have a = a* (k/K)~ 2 . 
Thus, following the above expressions, the temperature of the universe at the time when the 
mode k leaves the Hubble radius can be determined as 



T(r k ) ~ T. ( tt ) , (55) 



where k* and T* are associated with the initial moment of thermal equilibrium as introduced 

H(r k ) 
a(r k ) 



in previous subsection. Similarly the Hubble parameter H{r k ) = ^ Tfc ? during that time can 



also be derived using Eq. ff40|) which reads as 

/ k \ ^ {u+1) 

H(r k ) ~ I- , (56) 

where H* again is associated with k*. Thus the above two equations yield 

k\ 2 / A; x 3 



T(r k )~T^-j , H(r k )~H^-j , (57) 

if the perturbation modes exit the Hubble radius before the universe is fully thermalized. 
Here, k* is a normalization scale which is associated with the initial moment of thermal 
equilibrium, T* is the corresponding temperature at that moment, and if* is the Hubble 
parameter corresponding to fc* as well. Thus the heat capacity is given by 

4vr T 3 / k\~ 3 

Cv(r k ) « — 9sUr k )^ H > ( 5 §) 



3 *^-"'H* \k, 

at the moment of Hubble crossing during initial radiation contraction. Then we get the 
amplitude of gravitational potential as 

(-1^)1/2 ft— 

= k\ (59) 
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where Eq. ( |49p was applied and the scale factor a* is associated with the moment of thermal 
equilibrium. Following Eq. (T4"3"j) . we see that the amplitudes of the constant £L mode and the 
growing 5_ mode during the initial radiation contraction are related to the metric fluctuation 

as 

d_ 0) ~ $ H , s4k)~(^y$ H . (60) 

Here "Huv is the maximal scale of the Hubble parameter, in the beginning of the bouncing 
phase, which is roughly of the same order of kyy introduced in Eq. (jHJ). Thus from the above 
equation and Eq. (144)) we easily find out that the amplitude of the D + mode in expanding 
phase is mainly contributed by the S_ mode in IR regime. Correspondingly, we derive the 
primordial power spectrum of metric perturbation in IR regime as 

fc3 in I 2 ~ 3 9srJiyyT^ ( k \ 



* 2tt 2 ' +l " lQwMffiHj \k , 
and the corresponding spectral index is given by 

d In P$ . , 

which implies that the power spectrum of metric perturbations seeded by thermal fluctua- 
tions during the onset of radiation dominated contracting phase is a blue spectrum. 



4.2.2 Generation of primordial power spectrum of thermal fluctuations in UV 
regime 

As we have pointed out before, the thermodynamics of the radiation fluid in the contracting 
phase would be greatly affected by the Lee- Wick resonances at very high temperatures 
(T > M) and the corresponding equation of state of the radiation fluid evolves from the 
conventional value of w r = 1/3 to a non-conventional value of w r ~ 1 when more and more 
Lee- Wick resonances thermalize with the increasing temperature of the universe. Therefore, 
the generation of primordial power spectrum of thermal fluctuations is dramatically changed 
in the UV regime as the background and the radiation fluid at these later stage of the 
evolution will evolve with state parameters as w ~ w r ~ 1. 

During the Hubble-crossing of the mode k the cosmic fluid would be dominated by 
Lee- Wick resonance infested radiation plasma and the energy density of that fluid can be 
expressed as given in Eq. ffT5]) (where we can neglect the contribution of the CDM component 
required for the non-singular bounce). Also, to let the Hubble radius of the universe contract 
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which allows the mode to become super-Hubble, the first term in Eq. (j!5p would be the 
dominant term in the total energy density (as discussed before). If during that time the 
mass of the heaviest thermalized Lee- Wick particle is M then using Eq. ffl5|) in Eq. f )50|) one 
gets an expression for the heat-capacity of the radiation fluid at the time of Hubble-crossing 
of the mode k as 

C v (r k ) « V(r k )a(T k )M 2 T(r k ) , (63) 

where we have introduced 

/ \ g*N(T k ) 

= 12 , (64) 

which is associated with the number of internal degrees of freedom of the thermal system at 
the time r k . As the energy density and pressure of the radiation fluid is dominated by the 
a term (as given in Eq. ( TT5l) ) during the Hubble-crossing of the modes the equation of state 
of that plasma would be w r (Tn) ~ 1. In that case, following Eq. ( )63l) . the heat capacity of 
the universe would be 

Cv{r k ) - f «(r fe )M 2 |^ (£) \ (65) 

where we have used Eq. (1551) and Eq. fl56l) . Thus, if we also consider in this scenario that the 
background cosmology of the universe is also evolving due to the dominance of the Lee- Wick 
infested radiation field and consider w = w r ~ 1 during the Hubble-exit of the mode k, then 
we get the amplitude of the metric fluctuation corresponding to mode k using Eq. ( 1491) as 

= ^ ^ . (66) 

where a* is the scale factor corresponding to k* through fc* = a*H*. Following Eq. (166 p . 
we see that the amplitudes of the constant D_ mode and the growing £*_ mode during the 
contracting phase are related to the metric fluctuation as 

S^(k)^(^^j (67) 

where v = 2 as w = 1 in this particular case. Here "Huv is the maximal scale of the Hubble 
parameter, in the beginning of the bouncing phase, which is roughly of the same order of 
/cuv introduced in Eq. fj44l) . Thus from the above equation and Eq. (I44p we see that both 
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the modes D_ and S_ contribute equally to the amplitude of the D + mode which is the 
dominant mode of the expanding phase (as the other S + mode is the decaying one and 
its amplitude falls of with timejj. Hence, we obtain the primordial power spectrum of the 
expanding phase of the universe as 

p * 3 ln ,2 9.nM 2 t; fk\ 3 

F *=2^ |i,+I ^SS^ftUJ ' (68) 

and the corresponding spectral index is given by 

d\nPq> 

"*" 1 + 7taF = 4 - (69) 

which implies that the power spectrum of metric perturbations seeded by thermal fluctu- 
ations during the Lee- Wick resonance dominated contracting phase is highly blue in UV 
regime. 

4.2.3 Stability analysis and a rough estimate of the mass bounds of the Lee- 
Wick partners of chiral fermions and massless gauge bosons 

As the primordial power spectrum of metric perturbation generated in Lee Wick radiation 
phase is deeply blue, its amplitude would become secondary in infrared limit which cor- 
responds to large length scale. However, one needs to be aware of the potential concern 
that the power spectrum would become too large in ultraviolet regime where k takes a large 
value. Fortunately, for all nonsingular bounce models, there is a natural ultraviolet cutoff on 
k modes due to the existence of the bouncing scale. This is because, if a perturbation mode 
has not yet evolved into super-Hubble scale before the bouncing phase, then it will never 
take place unless there is a period of inflation after the bounce. As a consequence, one can 
easily read that the maximal value of the power spectrum takes place right at the beginning 
of the bouncing phase. 

The universe is in thermal equilibrium, with standard particles, in the phase of traditional 
radiation domination and thus one can assume the benchmark value of the Hubble parameter 
to be, H* ~ | j=^Z where g is the number of internal degree of freedom of traditional 
radiation. Using this equation and the relation between H and k from Eq. (|56p . one can 
calculate the maximal value of the primordial power spectrum as: 

pmax „ 27g^H^ 
* ~ 64vn ? 3 / 2 M p Ti ! ' { } 



2 We would like to refer to Refs. p21 [22l [26l [27l [29] for extensive analysis on microscopic description of 
nonsingular bouncing phase in a wide class of linear bounce models. 



22 



In the above equation i^uv is roughly similar to fcuv whose value is set by the microphysics 
of the bounce. In our notation it is trivial to show that 



where P$(k) = A; 3 |$fc| 2 /27r 2 is the power spectrum of metric perturbation $. Note that, this 
UV cutoff can be selected as the Planck scale /cuv ~ aM p in inflationary cosmology, and due 
to the quasi exponential expansion of the background universe, it leads to the well-known 
Trans-Planckian problem for inflationary perturbation [301 [31] ( see a l so [22] )• This issue does 
not happen in this model of bouncing cosmology as there is another natural cutoff, fcuv as 
shown in Fig. [2] which is much lower than Planck scale, i.e., the bounce scale (the maximal 
absolute value that the Hubble parameter can reach throughout the whole evolution). If 
P$(k) oc k n (where n = 3 in our bounce model) then we obtain, (<3>(x) 2 ) ~ ~ — — . As 
$(x) is supposed to be metric perturbation we expect ^($(x) 2 ) to be smaller than unity. 
Consequently the above equation implies the maximal value of the power spectrum to be 
less than unity as well. 

If we take g ~ 107 which is the number of internal degree of freedom of standard model 
of particle physics, then we can roughly obtain an upper bound on a combination of the Lee 
Wick mass and the bounce scale as 



by requiring P| iax ~ 1- We further take g*N ~ 100 and assume T* ~ 100 GeV, and then get 
M\Hw\ < 10 14 GeV 2 . However, since the bounce scale must be higher than the mass scale 
of Lee Wick partners |-£/uv| > M, it roughly implies that M < 10 7 GeV. 

This theoretical constraint restricts the masses of the Lee- Wick partners in the expanding 
or in the contracting phase of the universe to be less than 10 7 GeV as only the terms 
proportional to a is taken into account in the heat capacity. The constraint does not work 
in the bouncing phase of the universe (in the time between and r^") where much heavy 
Lee- Wick partners can get thermalized. In the present models where a dark component 
is used we have calculated the maximum mass of the Lee- Wick partner to be of the order 
of 10 15 GeV, but this does not contradict the analysis given above as the heavily massive 
Lee- Wick particle is only thermalized at the bounce point, near tb- Predominantly in the 
bouncing phase the Lee- Wick partners of the chiral fermions and mass less gauge bosons will 
contribute and consequently we can infer that most of the Lee- Wick partners of the chiral 
fermions and mass less gauge bosons will have their mass in the range 10 7 GeV — 10 15 GeV 
if the temperature near the bounce is 10 16 GeV. 




(71) 




(72) 
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4.2.4 Numerical analysis of the background and the perturbation in the bounc- 
ing phase 

Previously the problem of mode matching near the bounce was done analytically for modes 
in the IR or UV regime. In the analytic way one matches the modes which were super 
Hubble before tg and after tg. The super Hubble modes briefly re-enter the causal horizon 



during the bouncing phase (in the time interval between tg and t B ) 
inside the bouncing phase we use numerical techniques. 



To see what happens 
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Figure 3: Figure showing the behavior of the scale factor a(t) in the bouncing phase. 



To ensure that the model under consideration indeed yields a bouncing phase of the 
Universe we analyze the behavior of the background plasma on which the perturbations live. 
It is seen that with some benchmark values of the theory the model gives rise to predictable 
bouncing behavior. In the numerical code we set the scale factor a(t) = 1 at the bouncing 
time t = ts = 0. We also set the temperature of bounce, Tb, and the mass of the heaviest 
Lee- Wick partner's mass during the bouncing phase, M, as 



M = 10" 5 M p , T B = 10~ 4 M 



v ■ 



where M p is the Planck mass which is set to one, 

time unit is set to be the Planck time. The other parameters appearing in Eq. (E 

as 



(73) 



M p = 1. In this numerical analysis the 

are set 



g*N 



120 



9f 



40 



n 



10. 



(74) 
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Figure 4: Figure showing the behavior of the Hubble parameter H(t) in the bouncing phase. 

It is to be noted that as we treat the CDM component to be non-relativistic throughout the 
evolution of the universe, the energy density of this CDM component evolves as 



where the constant ps is found from the condition that at the bouncing time t = ts one 
must have Pr^b) + Poits) = 0. 

With these typical parameter values it is seen that the background properties of our model 
encoded in a(t), Hit), T(t), pit) and p(t) show a continuous non-singular bounce. From Fig. [3] 
it is clear that the scale factor a(t) has a minimum at t — 0. Before t = it is seen that a(t) 
is decreasing specifying a contracting universe which smoothly transforms into an expanding 
universe after the bounce. In Fig. H] the variation of the Hubble parameter H(t) across 
the bounce is depicted. The variation of temperature of the background near the bouncing 
point is shown in Fig. [5j If one compares the variation of energy density with time, as shown 
in Fig. O with the variation in temperature near the bouncing point one notices that the 
temperature of the system is the maximum at the bounce point when the total energy density 
is the minimum. This point shows the nature of Lee- Wick particle dominated cosmology 
where the Lee- Wick partners contribute to the energy density of the cosmic plasma with 
a negative contribution. In Fig. [7] the pressure of the background plasma, p(t), is plotted 
against time. The figure clearly shows that pressure is negative in the bouncing region, which 
is a prerequisite of the bouncing model. 
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Figure 5: Figure showing the behavior of the background temperature T(t) in the bouncing 
phase. 




Figure 6: Figure showing the behavior of the total energy density p(t) in the bouncing phase. 



The metric perturbation amplitude |$fc| for three values of the wave number are plotted 
in Fig. |HJ Numerical analysis shows that the perturbations remain almost constant across the 
bouncing point. The perturbations smoothly evolves through the bouncing point and that 
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Figure 7: Figure showing the behavior of pressure p(t) in the bouncing phase. 
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Figure 8: Figure showing the behavior of the perturbations across the bouncing point. In 
the above figure we have numerically evolved the perturbations along the bounce point for 
the following values of the wave number: from top the wave numbers are k = 8 x 10~ 7 , 
8 x 10~ 9 , 8 x 10" 11 respectively. The y— axis is plotted in logarithmic scale. 



27 



allows one to match the mode functions in the expanding phase after the bounce. While 
comparing these different Fourier modes, one finds that the perturbation mode with the 
highest value of k has the largest amplitude. Thus, it shows explicitly a blue spectrum is 
achieved in our model based on thermal initial condition. Therefore, the numerical result is 
in agreement with the analytic analysis performed previously. 

4.3 Thermal fluctuations in the case which does not include CDM 

Here we will discuss the scenario where a bounce can be achieved in the contracting phase 
without invoking a CDM component in the cosmic plasma. The main feature of this model 
is that the bounce would occur before all the heaviest Lee- Wick resonances could thermalize. 
So, some of the standard model particles will contribute to the radiation plasma without 
their Lee- Wick resonances (like the standard radiation phase) near the bounce. In this case 
the dominant part of energy density before or after t% is assumed to be coming from 
the standard radiation part, namely from cT 4 /4, in Eq. (I3"2"j) . The energy density starts 
to decrease (or increase) rapidly in the bouncing phase. Consequently in this scenario the 
dominant nature of the radiation plasma would be more like a conventional radiation fluid 
with w r ~ 1/3. Hence, for the modes which would leave the Hubble radius in the contracting 
phase, much before the bounce takes place, one can consider w = w r ~ 1/3. 

The analysis of the IR regime in the other case holds true here as well as the scenario 
of onset of radiation era are same in both the cases. We need to analyze the UV regime 
separately as now the background and radiation fluid evolve as w ~ w r ~ 1/3. Following 
the arguments given in the previous sub-sections, one can see that the heat capacity of the 
plasma during the Hubble-exit of mode k in this case would be 



where we have applied w = w r ~ 1/3 and consequently v = 3. In this case, by matching the 
formula (ITT]) and the general solution ( 143]) in contracting phase, the _D_ and S- modes can 
be expressed as 




(76) 



and the amplitude of the metric fluctuation with wavenumber k would become 




(77) 
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which show that the contribution of the D_ mode to the D + mode in the expanding phase 
is scale-free whereas the contribution of the S_ mode for the same is proportional to 1/k 
(following Eq. (jSJ)). Thus, if the S- mode (being the growing mode in the contracting 
phase) contributes dominantly to the D + mode in the expanding phase, then we get the 
power spectrum as follows, 

which is also a blue-tilted spectrum as before with a spectral index 

d In Pa, , . 

Note that, if the D + mode in the expanding phase is mainly inherited from the D_ mode in 
contracting phase for some reason, the corresponding power spectrum is proportional to k 3 . 

Thus from the above analysis we see that here the spectrum remains equally blue tilted 
in both IR and UV regime. In the other case, where a CDM component is added, the 
spectrum is more blue titled in the UV regime than that of IR. Therefore, it is not difficult 
to distinguish these two models observationally. 



5 Conclusion 

In this article an attempt has been made to apply the esoteric topic of Lee- Wick thermo- 
dynamics in the physics of the early universe undergoing a non-singular bounce. Bouncing 
cosmologies try to evade the Big Bang singularity which plagues most of the recent cos- 
mological models. Moreover the bouncing universe model presented in the article do not 
require inflation to solve the standard cosmological problems as an inflation-like scenario is 
built into the bouncing models. 

If one includes Lee- Wick partners of standard fields in the theory then it naturally leads 
to a radiation dominated phase where the partners of some of the standard fields contribute 
with a negative sign in the effective energy density of the universe. This negative radiation 
can lead to a bouncing phase of the universe when the total energy density reaches zero. 
There can be two options to produce the bounce. The first option includes a dark matter 
component which helps to produce a negative pressure near the bounce point[f|. A brief 
numerical analysis of the background evolution and the evolution of the perturbations across 
the bouncing point for this option is presented in the article. The other option does not require 

3 It helps to evade the null energy condition near the bounce point. 
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any dark matter component but it requires that the temperature of the universe near the 
bounce point to be less than the masses of some of the Lee- Wick partners of the standard 
fields. 

It has been shown that for the bounce mechanism to succeed one has to have the Lee- 
Wick partners of most of the chiral fermions and massless gauge bosons to be heavier than 
a mass scale. This scale depends on the bounce energy scale. This fact gives us a hint about 
the mass of the illusive Lee- Wick partners. 

In the bouncing universe scenarios presented in this article, the relevant perturbation 
modes leaves the Hubble radius while the universe is contracting as the Hubble radius starts 
to shrink. These modes lie outside the Hubble radius for a brief time. After the bounce is 
over and the Hubble radius starts to increase these modes starts to re-enter. In this way 
the bouncing universe scenarios can mimic the effects of inflationary universe models. The 
bouncing universe models naturally sets a wave number cut-off, in our case it is called fcuv- in 
the Lee- Wick thermal radiation induced bouncing models discussed here we have calculated 
the power spectrum of the metric perturbation $. In both the cases one gets a blue-tilted 
power spectrum. Therefore, the primordial power spectrum induced by thermal fluctuations 
in contracting phase is hardly to explain the CMB data since cosmological observations have 
proven a nearly scale-invariant power spectrum. However, this problem can be generally 
circumvented by introducing another light scalar field during matter contracting phase, which 
is known as the bounce curvaton mechanism [33l 131] . 

The theory with Lee- Wick partners are still in their infancy as they were neglected for 
a long time because there were other interesting scenarios in particle physics which could 
address the most pressing problems, as the hierarchy problem and other related problems, 
there. At the present moment Lee- Wick standard model and Lee- Wick thermodynamics 
related works are also trying to address the problems plaguing particle physics and cosmology. 
In this regard the present article discusses about an interesting property of the probable early 
universe utilizing the properties of the Lee- Wick theories. 
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A Transfer relation between the modes solutions in the 
contracting and expanding phases 

Here we will give a brief description of how to match the mode functions of the two phases, 
contracting and expanding, at the non-singular bouncing point by making use of the analysis 
of [2S]. In [22], the matching conditions of modes have been extensively derived in a situation 
where the stress-energy tensor undergoes a finite discontinuity at a phase transition. If this 
sudden change in the stress-energy tensor is due to a sudden change in the equation of state 
w, then the energy density remains constant on the hypersurface S of phase transition. In 
such a situation the two matching conditions in conformal Newtonian gauge (longitudinal 
gauge) read 



0. 



v k - Y (H> 







where v is the Mukhanov- Sasaki variable which can be written in terms of $ as 



3 (1 + w) 



h 2 -h 



7 ($' + n$), 



(82) 



Using r v = ^ the first matching condition yields 



D_ 



S—1~L— — SJHj 



(83) 



where we have used the solutions for <3>^r given in Eq. ( 1431) . On the other hand the second 
matching condition gives 



W_ - Hi 



U 1 14+ 

It i — 7L_ 



Hi + — D_ + 



ni + 



3 
k 2 



SJH. 



H'_ 
Hi 



Hi + - 
ui + 



k 2 

3 
~3 



Using the first relation given in Eq. ( 1831) in the above equation one gets 



2H' + - H 



D 4 



2Ui 



Hi - 
SJH. 



-h 2 _ e 

Hi + 3 



Hi 
1 



Hi H 



Hi 



nil n/2 nil nil 

rL_ — rL_ Tl i — Tl_ 



(84) 



55) 
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If the universe goes through a symmetric evolution around the bounce then the above equa- 
tion leads to 

D + = AD_ + Bk 2 S_ , (86) 
where A and B are constants of 0(1). 
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